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Abstract. Tropical cyclones are affected by a large number of climatic 
factors, which translates into complex patterns of occurrence. The variabil- 
ity of annual metrics of tropical-cyclone activity has been intensively stud- 
ied, in particular since the sudden activation of the North Atlantic in the mid 
1990's. We provide first a swift overview on previous work by diverse authors 
about these annual metrics for the North- Atlantic basin, where the natural 
variability of the phenomenon, the existence of trends, the drawbacks of the 
records, and the influence of global warming have been the subject of inter- 
esting debates. 

Next, we present an alternative approach that does not focus on seasonal 
features but on the characteristics of single events [Corral et al., Nature Phys. 
6, 693 (2010)]. It is argued that the individual-storm power dissipation in- 
dex (PDI) constitutes a natural way to describe each event, and further, that 
the PDI statistics yields a robust law for the occurrence of tropical cyclones 
in terms of a power law. In this context, methods of fitting these distribu- 
tions are discussed. 

As an important extension to this work we introduce a distribution func- 
tion that models the whole range of the PDI density (excluding incomplete- 
ness effects at the smallest values), the gamma distribution, consisting in a 
power-law with an exponential decay at the tail. The characteristic scale of 
this decay, represented by the cutoff parameter, provides very valuable in- 
formation on the finiteness size of the basin, via the largest values of the PDIs 
that the basin can sustain. We use the gamma fit to evaluate the influence 
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of sea surface temperature (SST) on the occurrence of extreme PDI values, 
for which we find an increase around 50 % in the values of these basin-wide 
events for a 0.49 °C SST average difference. 

Similar findings are observed for the effects of the positive phase of the At- 
lantic multidecadal oscillation and the number of hurricanes in a season on 
the PDI distribution. In the case of the El Nino Southern oscillation (ENSO), 
positive and negative values of the multivariate ENSO index do not have a 
significant effect on the PDI distribution; however, when only extreme val- 
ues of the index are used, it is found that the presence of El Nino decreases 
the PDI of the most extreme hurricanes. 
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1. Introduction 

Tropical cyclones are a rare phenomenon, with less than 100 occurrences per year world- 
wide (tropical depressions not counted) - just compare with 10 4 earthquakes with magni- 
tude larger than 4 per year. Despite this scarcity, the societal impact of these atmospheric 
systems is huge, as it is perceived year by year by the general public. An important issue 
for planning of tropical-cyclone damage mitigation is the knowledge of the temporal vari- 
ability of the phenomenon, and to try to establish which part of the variability is natural 
and which part can be due to global warming. 

When counting these meteorological monsters, a bit of terminology is useful [Emanuel, 
2005a]. Tropical cyclones is the generic name encompassing typhoons, hurricanes, tropical 
storms, and tropical depressions (when the context allows it, we will also use the vague 
term "storm"). Typhoons and hurricanes (and severe cyclonic storms) are "mature" 
tropical cyclones (the difference in name is only of geographical origin), with winds strong 
enough to achieve a category from 1 to 5 in the Saffir-Simpson scale [Schott et al, 2010]. 
In contrast, tropical storms have not so large wind speeds, and for tropical depressions 
these are even weaker. The thresholds separating these 3 classes are 34 and 64 knots for 
the so-called maximum sustained (1 min) surface (10 m) wind speed (1 kt = 0.5144 m/s), 
further thresholds define the Saffir-Simpson category. Then, speed is a way to define 
intensity, and intensity at a given instant determines the stage of the tropical cyclone 
(hurricane, tropical storm, etc.) and the category in case of hurricanes or typhoons. The 
stage corresponding to the maximum lifetime intensity determines how the storm is finally 
labeled. Generally, the databases in use only contain more or less complete information 
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on tropical storms and hurricane-like systems, and therefore, when talking about tropical 

cyclones, tropical depressions will be excluded in this text. Note that in some references 

tropical depressions are not even considered as tropical cyclones, they are just formative 

or decaying stages of tropical storms and hurricanes [Neumann et al, 1999]. 

Although the North Atlantic ocean (including the Gulf of Mexico and the Caribbean) 
comprises only about 10 % of the tropical cyclones in the world, this is the ocean basin 
which has been more extensively studied, with routine aircraft reconnaissance since 1944 
and an intensive scrutiny of ship log records for previous years (satellite imagery started to 
be available from the 1960's [Neumann et a/., 1999]). This has yield the longest and most 
reliable database for tropical cyclones, extending back in time, with obvious limitations, 
for more than one century and a half. 

In this work we analyze tropical cyclone activity in the North Atlantic. First, previous 
work on temporal variability is overviewed, from the point of view of direct observations, 
and the controversy about if a global-warming signal can be separated from the natural 
fluctuations is briefly mentioned. Fully in-depth reviews on this topic, including modeling, 
theory, and paleoclimatic studies, have been published by Shepherd and Knutson [2007], 
Knutson et al. [2010b] and Knutson et al. [2010a]. Second (Sec. 3), an approach based 
on the study of the features of individual tropical cyclones is presented, based on the 
calculation of their power dissipation index (PDI). Power-law fits for the PDI distribution 
are presented and discussed also, and a superior, more general fitting function (the gamma 
distribution) is introduced in Sec. 4, taking into account the exponential decay of the tail 
of the distribution. Sections 5 and 6 analyze the effect of tropical sea surface temperature 
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and diverse climatic indices (El Nino, NAO, AMO, number of tropical cyclones) on the 
PDI distribution. 

2. Variability of hurricane activity 

The large temporal variability of North-Atlantic tropical-cyclone frequency is clearly 
seen in the existing records, with a maximum of 28 occurrences in 2005 (including a 
subtropical storm) versus only 4 in 1983, for instance, see Fig. I. Excluding tropical 
storms, there was a maximum of 15 hurricanes also in 2005 and just 2 in 1982; and for 
major hurricanes, which are those with category 3 or above (96 kt threshold), there were 
8 in 1950 versus none in 1994 [Neumann et al, 1999; Gray et ai, 1992; Eisner and Kara, 
1999; Bell et al, 2006; Landsea, 2007a]. 

But high and low tropical-cyclone activity come in multidecadal clusters. Gray [1990] 
reported that the number of major hurricanes decreased to less than one half from the pe- 
riod 1947-1969 to 1970-1987 (3.3 per year versus 1.55), which was interpreted as a regional 
manifestation of a global-scale climate variation governed by thermoaline processes. If an 
hypothetical overestimation of speeds previous to 1970 is corrected, the decreasing trend 
still persists for these periods [Landsea, 1993]; nevertheless, see also Landsea [2005] with 
regard the necessity of this adjustment. 

Eisner et al. [2000] and Goldenberg et al. [2001] realized that in 1995 a new period of high 
activity seemed to have started in the North Atlantic. In the first case, the increase was 
related to the North Atlantic oscillation (NAO), whereas in the second it was associated 
to simultaneous increases in North Atlantic sea-surface temperatures (SST) and decreases 
in vertical wind shear, governed by variations of natural origin of the so-called Atlantic 
multidecadal mode (or oscillation, AMO, [Kerr, 2000]). In this way, the period from 1965 
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to 1994 witnessed no more than 3 major hurricanes per year, whereas in 1995 there were 
5 of such storms, and 6 in 1996. The return of high levels of activity was labeled as 
dramatic. 

Later, Webster et al. [2005] found that an increasing trend for the period 1970-2004 
was also significant for hurricanes of at least category 1, in clear contrast to the behavior 
of other basins worldwide (in which tropical SST also increased). Remarkably, these 
authors also extended the increase of major tropical cyclones to the other basins (excluding 
category 3) and noted that this was not inconsistent with climate simulations showing that 
a doubling of atmospheric C0 2 should lead to an increase in the frequency of the most 
intense tropical cyclones. However, subsequent analysis have claimed that the sign in the 
trend of the Northeastern Pacific is dependent on the time window selected and that the 
increase in other basins could be an artifact due to the lack of quality and homogeneity of 
the records. In any case, the increased activity of North-Atlantic hurricanes of categories 
4 and 5 was clear and robust [Klotzbach, 2006; Kossin et al, 2007]. 

However, Landsea [2007b] has shown that the analysis of tropical-cyclone variability 
using counts of events requires enormous caution. Before the advent of aircraft reconnais- 
sance in the Atlantic, in 1944, detection of open-ocean storms was dependent on chance 
encounters with ships. But even after the 1940's, aircrafts were covering essentially the 
west half of the basin, so systems developing on the east part were not always observed. 
The comparison between the number of landfalling storms (which is assumed reliable since 
1900) and the total number of storms gives a quite stable ratio of 59 % since 1966 (when 
satellite imagery started to operate), in contrast to a value around 75 % for previous 
years. In addition, new operational tools, as QuikSCAT, that have became available in 
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the turn of this century, have allowed the identification of one additional tropical cyclone 
per year (on average). This has led Landsea to estimate a deficit in the annual frequency 
of storms of 3.2 events per year up to 1965 and 1 per year from 1966 to 2002. (One can 
wonder why the correction to this problem is solved by adding a constant term, rather 
than multiplying by some factor.) 

Obviously, storm counts do not provide a complete characterization of tropical-cyclone 
activity, and other indicators are necessary. (Here we will use the term activity in a 
broad sense, and not as a synonymous of frequency or abundance.) In this way, total 
storm days, defined by Gray et al. [1992] as the total time all storms in a year spend 
in a given storm stage, was found by Gray [1990] to quadruple between low and high 
activity periods for the stage of major hurricanes (a reduction to 2.1 days per year from 
a previous average of 8.5, for the periods mentioned above). Note that this metric is 
different to the total duration of major hurricanes (which is the sum of all their durations, 
in any stage). For the total duration of North Atlantic hurricanes (categories 1 to 5, and 
including tropical-storm stages), Webster et al. [2005] reported a significant trend from 
1970 to 2004. 

Another quantity that has been employed is the net tropical cyclone activity (NTC), 
explained in the supplementary information of Goldenberg et al. [2001] as an average 
between storm counts and total storm days, considering tropical storms, non-major hur- 
ricanes, and major hurricanes, but giving more weight to latter and less weight to the 
former storms. In this way, information on frequency, duration, and intensity is combined 
into a single number, resulting that the NTC between 1995-2000 is twice the value of 
1971-1994, which in turn was a factor 1.6 smaller than for 1944-1970. 
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A significant step forward was taken by Bell et al. [2000] , who defined the accumulated 
cyclone energy (ACE) by summing, for all tropical cyclones in a season, the squares of the 
6-hour maximum sustained wind speed for all records in which the systems were above 
the tropical storm limit (i.e., speed > 34 kt). The ACE is a more natural way to combine 
frequency, intensity and duration than the NTC, but it should not be interpreted as a 
kinetic energy (rather, it could be the time integral of something akin to kinetic energy, 
as its name denotes). It is worth mentioning that previously, Gray et al. [1992] had 
introduced a related quantity, the hurricane destruction potential (HDP), with the only 
difference that tropical-storm stages were not taken into account. 

Updating the results of Bell et al. for the ACE as an annual indicator, Trenberth [2005] 
noticed that between 1995 and 2004 the North Atlantic activity had been "above normal", 
with the only exception of the El Nino years of 1997 and 2002, and that this increase was 
in parallel to the fact that that decade was the one with the highest SST on record in 
the tropical North Atlantic (by more than 0.1 °C) and also to the increased water- vapor 
content over the global oceans. (The results for ACE were confirmed later by Klotzbach 
[2006].) Although many uncertainties remain on how this and other environmental changes 
can affect the self-organization process of tropical-cyclone formation, once formation has 
taken place it seems clear that these conditions provide more energy to the storm, which 
suggests more intense winds and heavier rainfalls under a global warming scenario. 

A similar analysis was performed by Emanuel [2005b], who introduced what he called 
power dissipation index (PDI), but that we will call annual or accumulated PDI (APDI), 
defined as the ACE with the main difference that the square of the speeds is replaced 
by their cube (and also, the time summation was performed for the whole life of the 
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storm and not only for tropical-storm and hurricane stages). It is notable that the APDI 

constitutes a "proxy" estimation of the kinetic energy dissipated by all tropical cyclones 

in one season (not of their power). Emanuel showed how a well-known formula for the 

dissipation, used previously by Bister and Emanuel [1998], could be converted into the 

APDI formula under reasonable assumptions. Indeed, the dissipated power is given by 

the integral over surface of the cube of the velocity field (multiplied by the air density 

and the drag coefficient, which can be assumed as constants). First, accepting a similar 

shape for all storms (the same functional form for the velocity field, scaled by the radius 

of the storm and the maximum-in-space instantaneous speed) this power is proportional 

to the square of the radius multiplied by the cube of the maximum speed (with the same 

proportionality constant for all storms). Second, noticing the weak correlation between 

storm dimensions and speed, one could assign a common radius to all storms (in all stages) 

and get only random errors in the estimation of the power. As the energy dissipated by 

a tropical cyclone is given by the time integral (for all its life) of its power, replacing the 

integral by a discrete summation (as the records come in discrete intervals) one gets that 

the dissipated energy is roughly proportional then to the sum of the cube of the maximum 

sustained wind speed, which is (in contrast to the usage of Emanuel), what we call PDI. 

Further summation for all storms provides the APDI, for which it is expected that the 

errors coming from the assignment of a common radius will more or less compensate. See 

Corral [2010] for a non-verbal version of this derivation, and the appendix here for an 

estimation of the proportionally constant linking PDI and dissipated energy. 

In equations, 

n 

APDI = J2 PDh 

i=i 
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where % counts the n tropical cyclones in the season under consideration. For a single 
storm, its PDI is 

Mt 

where t labels discrete time, v t is the maximum sustained surface wind speed at t, and we 
introduce a time interval At that is constant, in principle, and equal to 6 hours, just to 
make the PDI independent on changes on At (nevertheless, it is an open question how a 
better time resolution can influence the stability of the PDI value, this will depend on the 
smoothness of the time-evolution of the speed). Under these assumptions, the SI units of 
the PDI are m 3 /s 2 (for instance, hurricane Katrina [2005] yields PDI= 6.5 ■ 10 10 m 3 /s 2 ). 
Note that for tropical cyclone i, 

PDIi = (vf)iTi, 

where {vf)i is the average of the cube of the (6-hour) maximum sustained wind speed over 
storm i, and T { is its duration. Also, for year y, 

APDI y = n y (PDI) y = n y (T) y (v*) y 

where n y is the number of tropical cyclones in that year and is the average of the 
considered quantity for the same year (the average is performed in a different way for 
the duration, of which there are n y data, than for the speed, of which there are J2iTi/At 
records) . 

The previous equations show that the PDI of a tropical cyclone is the product of its 
duration and its average intensity, if one redefines intensity as the cube of the maximum 
sustained surface wind speed, whereas the APDI is the product of frequency, average 
duration, and annually averaged (6-hour) intensity. Both the PDI and the APDI turn 
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out to be very natural and convenient ways to estimate tropical-cyclone activity, as they 
are a rough estimation of individual and total dissipated energy, respectively. Further, 
they are more robust than other measures of storms, as duration or track lengths, due to 
the fact that for these latter, the value is highly influenced by the definition of when a 
tropical cyclone starts, whereas for the PDI and APDI such values have little influence 
on the final result, as they are weighted by the smallest values of vf. 

What Emanuel found was a clear correlation between "unprecedented" increases in 
tropical SST and APDI, both for the North Atlantic (more than doubled in the last 30 
years) and Northwestern Pacific (75 % increase), when a 1-2-1 filter was applied twice 
to both signals (i.e., a 3- year running average where the central point has a weight that 
is equal to that of both neighbors). As the upswing in SST has been generally ascribed 
to global warming, this author argued that the increase in APDI could be partially of 
anthropogenic origin. 

However, Landsea [2005] noticed that, after applying the smoothing procedure, Emanuel 
failed to drop the ending points, which were not averaged and were substantially larger 
than the smoothed series for the North Atlantic, creating the impression of a dramatic 
increase in PDI values. Together with Gray [2006], Landsea has also criticized that 
Emanuel reduced the values of the wind speeds before the 1970's by an excessive amount 
(up to more than 20 knots for the largest values), in order to correct an overestimation in 
those values. In fact, recent research has shown that surface winds in major hurricanes are 
stronger than previously assumed, and therefore, probably no correction at all is necessary 
[Landsea, 2005]. The North Atlantic PDI series, with no correction, is plotted in Fig. 2. 
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Other problems with the accuracy of the records, particularly for the Eastern Hemisphere 
are pointed out by Landsea et al. [2006]. 

Many more articles have been devoted to these complex affairs, which unfortunately 
cannot be abstracted in this brief overview. Before ending this section, let us just mention 
the work of Chan [2006], Eisner et al. [2006], Emanuel [2007], Wu et al. [2008], Swanson 
[2008], Eisner et al. [2008], Aberson [2009], and Landsea et al. [2010], for its relation with 
the issues discussed here. 

3. PDI distribution and power-law fits 

All the studies summarized so far, counting numbers of storms, total number of storm 
days, calculating NTC, ACE, or APDI, have paid attention to overall measures of annual 
tropical-cyclone activity, that is, they were trying to answer the important question of 
comparing the characteristics of different seasons, or longer periods. But how different are 
the tropical cyclones from one season to another? i.e., which are the features of individual 
tropical cyclones in a given phase of activity, or under the influence of a certain value of 
a climatic indicator? 

The study of Webster et al. [2005] contains a first attempt in this direction, comparing 
the ratio between tropical cyclones in a given Saffir-Simpson category and tropical cyclones 
in all categories (1 to 5). For global aggregated data it was found that the proportion of 
category 4+5 events increased from less than 20 % in the early seventies to about 35 % after 
2000, with a corresponding decrease in the proportion of category 1 storms, from more 
than 40 % to 30 %. This is an indication that more major hurricanes (or typhoons) are 
present, in comparison, but not that the storms are becoming more intense, individually 
(i.e., there are more major ones in proportion, but their intensity is not necessarily record- 
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breaking). It is an open question in what part these results could be an artifact due to 
the incompleteness of the data for the earlier years [Gray, 2006; Klotzbach, 2006; Kossin 
et al, 2007]. In the North Atlantic, where the records are the best, these changes were 
much more modest, from a 20 % of category 4+5 in 1975-1989 to 25 % in 1990-2004. In 
fact, what Webster and coworkers were really doing was the calculation of the probability 
of hurricanes in a given category, and how this probability distribution changes from one 
period to another. 

However, in order understand changes in the behavior of individual tropical cyclones 
we need to know, first of all, which are the general properties of tropical cyclones, in 
other words, which is their "unperturbed" nature, and only after that we could study the 
influences of year-to-year variability on them. 

This has been attempted by one of the authors and collaborators [Corral et al, 2010], 
analyzing the statistics of the (individual-storm) PDI for long periods of time. As we have 
mentioned, the PDI is an approximation to dissipated energy and is therefore perhaps the 
most fundamental characteristic of a tropical cyclone. 

In this way, the PDI for each tropical cyclone in the period and basin considered was 
calculated, and the resulting probability density D(PDI) was obtained, following its 
definition, as 



where there is a certain freedom in the choice of the width of the interval (or bin) [Her- 
garten, 2002] and Prob denotes probability, estimated by relative frequency of occurrence. 
Of course, normalization holds, / °° D(PDI)dPDI = 1. For more concrete details on the 
estimation of D(PDI), see the supplementary information of Corral et al. [2010]. 



D(PDI) 



Prob [value is in an interval around PDI] 



width of the interval 
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The results for the 494 tropical cyclones of the North Atlantic during the period 1966- 
2009 are shown in Fig. 3, in double logarithmic scale. The data are the best tracks from 
NOAA's National Hurricane Center [Jarvinen et al, 1988; NHC, 2011]; as we have seen, 
1966 marks a year from which the quality and homogeneity of the records is reasonable 
(disregarding the possible overestimation of wind speeds previous to 1970). A straight line 
regime is apparent in Fig. 3(a), signaling the existence of a possible power-law behavior 
over a certain range of PDI values, 



where oc indicates proportionality, r is the exponent and a and b are the truncation 
parameters. 

Indeed, statistical tests support the power-law hypothesis. The original reference [Corral 
et al, 2010] used a generalization of the fitting and testing procedure proposed by Clauset 
et al. [2009]. The key point is to find the optimum values of a and b, for which the 
following steps are followed: 

• Fix arbitrary initial values of the truncation parameters a and b, with the only 
condition that b/a > 10, for example. 

• Find the maximum likelihood estimation of the exponent r, by means of the maxi- 
mization of 



D(PDI) oc 



1 



, for a < PDI < b, 



PDP 



£(r;a,b,G) = In 



t-1 



— t In In a, 

a 



1 - (a/b) 
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or, alternatively, the solution of 



r — 



1 




with G the geometric mean of the values comprised in the interval [a, b). 
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• Calculate the Kolmogorov-Smirnov (KS) distance between the empirical PDI dis- 
tribution (defined only for the range a < PDI < b) and its power-law fit [Press et al, 
1992]. 

• Look for and select the values of a and b that minimize the Kolmogorov-Smirnov 
distance. Select also the corresponding r exponent. 

Once the optimal values of a and b, and from here the value of r, are found, a p— value 
quantifying the goodness of the fit can be calculated just by simulating synthetic data 
sets, as close as possible to the empirical data (with the selected power-law exponent r 
between a and b and resampling the empirical distribution outside this range). The pre- 
vious 4 steps are applied then to each of these synthetic data sets (exactly in the same 
way as for the empirical data, in order to avoid biases), which yields a series of minimized 
Kolmogorov-Smirnov distances, from which one can obtain the p— value, which is the sur- 
vivor probability function of the synthetic KS distances evaluated at the empirical value. 
To be more precise, p is defined as the probability that, for true power-law distributed 
data, the KS distance is above the empirical KS value. In this way, very small p— values 
are very unlikely under the null hypothesis, and one must conclude that the data are not 
power-law distributed. For technical details about the whole procedure, see the supple- 
mentary information of Corral et al. [2010]. The introduction of the (arbitrary) condition 
b/a > 10 is necessary because in the case of a double truncation (from above and from be- 
low) several power-law regimes can coexist in the data (in contrast to the original Clauset 
et al.'s case, b — >■ oo), and one needs a criterion to select the most appropriate one, in our 
case that the range is high enough. 
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The results obtained from the application of this method to the North-Atlantic 1966- 
2009 PDIs appear on table 1, and they show that there is no reason to reject the power- 
law hypothesis for a certain set of values of a and b. Among the several outcomes of the 
method, depending on the conditions imposed, a good choice is the one with the values 
of a and b around 2 • 10 9 and 6 ■ 10 10 m 3 /s 2 , with a resulting exponent r = 1.16 and a 
p— value about 40 %, see Fig. 3. In general, one needs to balance the tendency to select 
low b/a ratios, which yield high p— values, with the tendency to larger ratios and lower 
p. 

In Corral et al. [2010] it was checked that the behavior of the PDI distribution is 
very robust, showing very little variation for different time periods, or if extratropical, 
subtropical, wave and low stages are not taken into account, or if U.S. landfalling tropical 
cyclones are removed from the record, or speeds previous to 1970 are reduced by an amount 
of 4 m/s [Landsea, 1993]. In addition, the power-law behavior was also found in other 
basins: the Northeastern Pacific, the Northwest Pacific and the Southern Hemisphere as 
a whole (the North Indian ocean was not analyzed due to the low number of events there 
and the lack of a complete satellite monitoring until recently [Kossin et al, 2007]). The 
resulting exponents were r = 1.175, 0.96, and 1.11, respectively, ±0.05, in the worst case. 

Further, it has been pointed out that the power-law PDI distribution could be a re- 
flection of the criticality of hurricane occurrence [Corral, 2010], in the same way as for 
earthquakes, forest fires, rainfall, sandpiles, etc. [Bah, 1996; Turcotte, 1997; Jensen, 1998; 
Sornette, 2004; Christensen and Moloney, 2005]. In fact, tropical-cyclone criticality could 
be only a part of the criticality of atmospheric convection [Peters and Neelin, 2006, 2009]. 
The parallelisms between tropical-cyclone occurrence and self-organized criticality are dis- 
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cussed in depth in Corral [2010]. The implications of this finding for the predictability 
of the phenomenon would be rather negative, in the sense that this kind of processes are 
characterized by an inherent unpredictability, although the differences with deterministic 
chaos are worth of being investigated in the future. 

Nevertheless, the previous statistical method was not found to be fully satisfactory 
when applied to a different data set [Corral et al, 2011]. The problem lies in the selection 
criterion for a and b. First, it can provide non optimal values of them, which would lead 
to the rejection of the null hypothesis in cases in which is true. Second, it has a strong 
tendency to underestimate the upper limit b (if r > 1). 

Another method, with a distinct selection procedure for a and b, was introduced by 
Peters et al. [2010]. The steps are: 

• Fix arbitrary initial values of the truncation parameters a and b, with no restrictions. 

• Find the maximum likelihood estimation of the exponent r (by means of the same 
formula as for the previous method). 

• Calculate the Kolmogorov-Smirnov distance between the empirical PDI distribution 
(defined only for the range a < PDI < b) and its power- law fit. 

• Calculate the p— value of the fit, which we denote by q, by simulations of synthetic 
power-law data sets, with exponent r (the simulation of the distribution outside the 
interval [a, b) is not necessary in this case). 

• Look for and select the values of a and b which maximize the logarithmic range, b/a, 
under the condition, let us say, q > 20%. 

In the original reference by Peters et al. the number of data in between a and b, 
N a b, was the quantity which was maximized, rather than the log-range b/a. As in our 
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case the power-law exponent is very close to one (for which there is the same number of 
data for each order of magnitude), there should be no relevant difference between both 
procedures; nevertheless, the we think the present choice is more appropriate for tails 
with larger exponents. Of course, setting the threshold q equal to 20 % is arbitrary, and 
we will have to check the effect of changing this value. 

One has to take into account that the q value calculated in this method is not the true 
p— value of our fitting procedure. The former is calculated for fixed, or known, a and b 
values, whereas the latter should take into account that both parameters are optimized. 
Comparing with the previous method it seems that the change between both p— values 
can be around a factor 2 or 3 (so caution is necessary at this point). The results of this 
method for the same North Atlantic data used in table 1 are displayed now in table 2. 
Indeed, a comparison between both tables show that the results are consistent, if one 
takes into account the difference between p and q. In summary, the existence of a power 
law over a range larger than one decade is well supported by the statistical tests. 

4. Gamma distribution of PDI 

So far, we have established that a "significant" range of the PDI probability density can 
be described as a power law. Deviation at small values (PDI < a) is justified through the 
fact that the best-track records are incomplete - this is obvious, as tropical depressions are 
deliberately excluded from the database. Inclusion of tropical depressions in the analysis 
of the Northwestern Pacific (which are easily available, in contrast to tropical depressions 
at the NHC best-track records) enlarges somewhat the power-law range, although the 
coverage of tropical depressions in that basin is far from exhaustive [Corral et ai, 2010]. 
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The rapid decrease of the PDI density above the value of b has been explained as a 
finite size effect: tropical-cyclones cannot become larger (in terms of PDI) because they 
are limited by the finiteness of the basin; either they reach extratropical regions (cold 
conditions) or they reach land, in which cases they are dissipated as they are deprived 
from their warm-water energy source. In this way, it has been shown how the tracks of 
most events with PDI > b are affected by the boundaries of the basin, and a finite-size 
scaling analysis showed how a reduced area (limiting it in terms of longitude) moved the 
cutoff to smaller values, with the cutoff defined roughly as the point in which the PDI 
probability density clearly departs from the power law [Corral et al, 2010]. 

Interestingly, it has also been found in the previously mentioned reference that the 
season-averaged tropical SST of the North Atlantic has an influence on the PDI values 
similar to a finite size effect. Years with high SST lead to a larger cutoff, just the opposite 
as years with low SST, but keeping nearly the same value of the power-law exponent (the 
same has been found also for the Northeastern Pacific). So, the power-law is a robust 
feature of the PDI distribution, but it is not telling us anything about the influence of 
other factors (finiteness of the basins and SST); rather, it is the cutoff, which is not 
properly defined yet, which carries this information. 

In order to overcome this deficit, a function modeling the PDI distribution that includes 
the tail seems appropriate. Experience with finite size effects in critical phenomena (taking 
place in continuous phase transitions or in avalanche-evolving systems) suggests a simple 
exponential factor for the tail, and therefore a gamma-like distribution for the whole 
domain (excepting the incompleteness for small values), 




DRAFT 



October 5, 2012, 1:29am 



DRAFT 



CORRAL AND TURIEL: VARIABILITY OF N. ATLANTIC HURRICANES 



X - 21 



where x plays the role of the PDI (only for aesthetic reasons), c represents now the cutoff, 
which is a scale parameter, and j3 is a shape parameter, with the power-law exponent 
equal to 1 + j3. If a — 0, j3 has to be negative, in order to avoid the divergence of the 
integral of D(x); however, for a > 0, there is no restriction on {3. The (complement of the 
unrescaled) incomplete gamma function is 



[Abramowitz and Stegun, 1965], with 7 > for z a = but unrestricted for z Q > 0. 
Unfortunately, some of the numerical routines we will use are not defined for 7 < 0, so, 
we will need to transform 



(integrating by parts), where, for j3 < 1 the incomplete gamma function on the right hand 
side can be computed without any problem. An alternative approach for the evaluation 
of the cutoffs would have been to compute the moment ratios used by Peters et al. [2010]. 

The reason to model D(PDI) by means of a gamma distribution is due to the fact 
that one can heuristically understand the finiteness of a system as introducing a kind of 
effective correlation length. Then, it is well known that, in the simplest cases, correlations 
enter into probability distributions under an exponential form [Zapperi et al, 1995]. 

In contrast to the method of fit and goodness-of-fit test explained above and to other 
previous work [Corral, 2009], in order to avoid difficulties, we perform a simple fitting 
procedure here, acting directly over the estimated probability density. The value of the 
minimum limit a is selected from information coming from the power-law fits (tables 1 and 
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2) complemented by visual inspection of the plots of the empirical probability density. A 
reasonable assumption is a = 2-10 9 m 3 /s 2 (note that the empirical D(PDI) is discretized). 

The fit of the shape and scale parameters, (3 and c, is performed using the nonlinear 
least-squares Marquardt-Levenberg algorithm implemented in gnuplot, applied to the 
logarithm of the dependent variable (i.e., it is the logarithm of the model density which 
is fit to the logarithm of empirical density); further, the parameter introduced in the 
algorithm is not c but d = lnc (i.e., we write c = exp(c')). Moreover, it is necessary 
to correct that the empirical density is normalized from to oo (i.e., not truncated, all 
events contribute to the normalization), whereas the model density goes from a to oo. 
In order that one can fit the other, a factor N a /N must multiply the fitting density (or, 
equivalently, divide the empirical one by the same factor), where N a is the number of 
data points (tropical cyclones) fulfilling x > a. Due to the discretization of the empirical 
density, in practice, a has to coincide with the lower limit of a bin of the density, or, 
alternatively, N a has to be redefined accordingly. 

The resulting fit for the 1966-2009 PDIs of the North Atlantic is displayed on Fig. 4, and 
the obtained parameters are r — 1 + /3 — 0.98±0.03 and c = (8.1 ±0.4) ■ 10 10 m 3 /s 2 , when 
the empirical density is estimated with 5 bins per order of magnitude (and the position of 
those bins is the one in the figure). In contrast with the method used for power-law fits, 
the results here depend on the width of the intervals of the empirical probability density 
(the number of bins or boxes per order of magnitude), and their position in the x— axis. 
We have found that, in general, the influence of these factors on the parameters is more 
or less within the error bars given above. 
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Note that the exponent r was larger than 1 for the (double truncated) power law, 
whereas now t — 1 + /3 turns out to be slightly below I. The difference is small, but, more 
importantly, both exponents represent different things, because we are fitting different 
functions. If the true behavior of D(PDI) were a gamma distribution (which is in practice 
impossible to know), it is expected that a pure power law (over a certain range) would 
yield a larger (steepest) exponent, due to the bending effect of the exponential factor. 

On the other hand, a gamma distribution with 1 + j3 < 1 cannot be extrapolated to the 
case in which the cutoff is infinite and cannot describe vanishing finite size effects, as it 
would not be normalizable. If we enforce 1 + /3 > 1 in the fit, we get 1 + /3 = 1.00±0.03 and 
c = (8.3±0.4)-l 10 m 3 /s 2 , which is not only a visually satisfactory fit but undistinguishable 
from the previous one. So, there is a tendency of the exponent to be smaller than one, 
but values slightly above one are equally acceptable. 

5. Sensitivity of the tail of the PDI distribution to SST 

We now investigate how the splitting of the PDI distribution into two parts, one for 
warm SST years and another for colder years, influences the shape and scale of the result- 
ing distributions. First, monthly SST with one degree spatial resolution are averaged for 
the tropical North Atlantic (90° to 20°W, 5° to 25°N) during the hurricane season (June 
to October). This yields a single SST number for each year or season [Webster et ai, 
2005], and in this way, values above the 1966-2009 mean are considered to denote warm 
years, and the opposite for cold years (notice that the comparison is not done relative to a 
different, longer period). The average SST for warm years turns out to be 0.49 °C larger 
than that of cold years. The SST data are those of the Hadley Centre, UK Meteorological 
Office [Rayner et al, 2003; BADC, 2010]. 
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The resulting PDI distributions are displayed on Fig. 5, which shows that the largest 

values of the PDI correspond to high SST years. If we compare the mean values of the 

PDI, the increase is around 40 %, i.e., (PDI) warm / (PDI) cold ~ 1.42 ± 0.19, which is 

significantly larger than 1 (at the 99 % confidence level, see table 3). 

The gamma fit also does a good job here, with an outcome, for high SST, 1 + j3 = 

0. 94 ± 0.06 and c = (9.6 ± 1) • 10 10 m 3 /s 2 , taking a = 2 • 10 9 m 3 /s 2 . For low SST, visual 
inspection of the plot of the density suggests that it is more appropriate to raise the value 
ofato3-10 9 m 3 /s 2 and then, 1 + /3 = 0.98±0.11 and c = (6.5±1) • 10 10 m 3 /s 2 . This yields 
a 50 ± 27% increase in the value of the cutoff, which is significantly larger than zero if one 
assumes normality for the c parameters. Note that we do not expect differences between 
the ratios of the cutoffs and the ratios of the means, due to the fact that the power-law 
exponent 1 + /3 turns out to be smaller than 1, and then the mean scales linearly with the 
cutoff, see the supplementary information of Corral et al. [2010]. 

If, in order to make the comparison of c in the same conditions between warm and cold 
years, we fix a = 3 • 10 9 m 3 /s 2 and 1 + f3 — 0.93 (this is the value obtained in this case for 
the whole dataset, with c = (7.8 ±0.4) • 10 10 m 3 /s 2 ), then we get c = (9.6 ±0.5) • 10 10 m 3 /s 2 
and (6.1±0.4)-10 10 m 3 /s 2 , for high and low SST years, respectively, which is not essentially 
different to the above case, but has smaller errors. This leads to c warm /c co id = 1.57 ± 0.13, 

1. e., more than a 50 % increase in the largest values of the PDI for an increase in the SST 
equal to 0.49 °C, which seems to indicate a high sensitivity of the extreme (cutoff) PDIs 
to the warming of the sea surface. The fits are displayed in Fig. 5. 

A further step is to separate the years not into two groups, warm and cold, but into 3, 
very warm, intermediate and very cold (in relative terms), and compare only the years with 
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"extreme" temperatures. Taking the thresholds as the mean SST plus half its standard 
deviation and the mean minus half the standard deviation we get a ratio of mean PDIs 
equal to 1.65 ± 0.25, for an increase in mean SST about 0.65°C, see table 3. Thus, the 
splitting of the distributions is more pronounced for "extreme" SST. 

6. Influence of El Nino, the NAO, the AMO, and the number of storms on 
the PDI 

We can proceed in the same way for other climatic variables, or indices. Considering 
the El Nino southern oscillation (ENSO), the multivariate ENSO index (MEI) separates 
warm phases of the tropical Eastern Pacific (MEI > 0) from cold phases (MEI < 0), with 
high values of the MEI associated to El Nino and low values to La Nina [Wolter and 
Timlin, 1998]. It has been proved that El Nino phenomenon is anti-correlated with hur- 
ricane activity in the North Atlantic, in such a way that the presence of El Nino partially 
suppresses this activity Gray [1984], which has been detected even in the record of hurri- 
cane economic losses in the U.S. [Pielke Jr. and Landsea, 1999]. Using data from Wolter 
[2010], and averaging for each season the months from May/June to September/October, 
we can distinguish between positive and negative phase seasons, and from here obtain the 
corresponding PDI distributions. 

We do the same for the NAO and AMO indices (North Atlantic oscillation and Atlantic 
multidecadal oscillation, respectively, averaging the indices from June to October), as well 
as for the annual number of tropical cyclones, utc, and the annual number of hurricanes, 
riH- In the latter case, only hurricanes are accounted in the average of the PDI. High 
values of MEI, NAO, and AMO refer to just positive values of the indices, whereas high 
values of utc an d uh denote values above the 1966-2009 mean (in the same way as for 
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the SST), and conversely for low values of the variables. Very high values correspond in 
all cases to values above the mean plus one half of the standard deviation (and conversely 
for very low values). 

The results can be seen in Figs. 6 and 7 and in table 3. Comparing PDI distributions 
for positive and negative values of the ENSO, AMO, and NAO indices, only the ones 
separated by the AMO show significant differences (at the 95 % level), with positive 
values of AMO triggering larger extreme PDI values. The lack of influence of ENSO in 
the PDI is in agreement with Corral et al. [2010]. 

But, if the PDI distributions are compared for high enough values of the indices, both 
the MEI and the AMO show a clear influence on the PDIs (not the NAO). As seen in Fig. 
6, the presence of El Nino (very high MEI), in addition to reduce the number of North 
Atlantic hurricanes (which is well known), decreases the value of their PDI. Of course, just 
the same happens for very low values of the AMO index, which decrease both hurricane 
numbers [Goldenberg et al, 2001] and the largest PDI values. La Nina or high values of 
the AMO have the opposite effect, increasing the most extreme PDIs. 

The influence of the number of tropical cyclones in the values of the PDI is not significant 
for values of n TC above or below its mean, but larger /smaller values of utc (mean ± 1/2 of 
standard deviation) are correlated with larger /smaller values of the PDIs, with an increase 
in (PDI) around 60 % when the mean utc goes from 7.35 events to 16.8. If we consider 
hurricane counts this effect is more pronounced, although care has to be taken with a 
kind of circular argument here: just by chance, the number of hurricanes can increase at 
the same time that the number of tropical storms decreases, then, it is likely that the 
mean PDI will be higher in this case (as the maximum sustained wind speed is correlated 
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with PDI), but no physical effect is behind this, only statistical fluctuations. Therefore, 
we study the influence of the number of hurricanes on the PDI distributions of hurricanes 
only, finding that years with more hurricanes also have larger PDI values. This seems to 
indicate that the conditions for genesis and survival of the storms are correlated; compare 
with the conclusions in Sec. 4.3 of Eisner and Kara [1999]. 

7. Summary and Conclusions 

We have illustrated how diverse metrics have been introduced by several authors in 
order to quantify tropical-cyclone activity, for instance: major- hurricane counts, total 
major- hurricane days, major-hurricane duration, NTC, etc... In general, these metrics 
show a clear increase in activity in the North Atlantic since 1995 -although it is difficult 
to associate a single cause to this phenomenon. Particularly interesting is the PDI, which 
was proposed by Emanuel [2005b] in order to estimate annual tropical-cyclone kinetic 
energy dissipation in different basins. 

In a previous work [Corral et al, 2010] we demonstrated how the use of the PDI to 
characterize the energy dissipation of individual tropical cyclones gave coherent and ro- 
bust results, leading to nearly-universal power-law distributions (i.e., with roughly the 
same exponents for different ocean basins). The outcome of this approach, in contrast 
to the findings obtained using other indices of activity, has important implications for 
the understanding of the physics of tropical cyclones, as it allows a connection with self- 
organized-critical systems [Corral, 2010]. This, in turn, opens interesting questions about 
the limits of predictability of tropical cyclones and, in a broader context, about the com- 
patibility between criticality and chaos in the atmosphere. 
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Due to the important implications that the emergence of power-law distributions has, 
we also discuss different ways of fitting them. The main problem is not to find the power- 
law exponent (which is the only parameter of the distribution) but to decide for which 
values of the variable the power-law regime holds. The procedure is based in the work 
of Clauset et al. [2009], although some variations seem necessary in order to improve the 
performance of the method. 

Whereas the power-law part of the PDI distribution does not change under different 
climatic conditions (as it should be in a critical system), the tail of the distribution 
does, containing therefore precious information about external influences on the system, 
in particular on the finiteness of the part of the basin able to sustain tropical-cyclone 
activity. Therefore, a probability distribution modeling not only the power-law regime 
but also the faster decay at the largest PDI values seems very appealing. We propose the 
use of the gamma distribution, combining a power law with an exponential decay. 

In this way, it is shown that the cutoff parameter modeling the exponential decay 
increases with annually averaged SST, which means that high SST has an effect on the 
PDI that is analogous to expanding the effective size of the part of the ocean over which 
tropical cyclones develop (and the opposite for low values of the annual SST). The same 
effect is found for the AMO index, with high AMO values leading to larger extreme PDI 
values. However the El Nino phenomenon has the opposite effect, being the presence of 
La Nina which triggers the largest (in dissipated energy) North Atlantic hurricanes. The 
number of hurricanes is also found to be positively correlated with their PDI. In contrast, 
for the NAO index we do not find any significant correlation. 
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In conclusion, the characterization of PDI probability distribution reveals as an inter- 
esting tool not only to learn about the fundamental nature of tropical cyclones but also 
to evaluate the effect of different climatic indices on the energy dissipated by them. 

Appendix: Relation between PDI and Dissipated Energy 

The power dissipated at any time by a tropical cyclone can be estimated by means of 
the formula [Bister and Emanuel, 1998], 



where t is time, r is the spatial coordinate over the Earth's surface, p is the air density, 
Cd is the drag coefficient, and v is the modulus of the wind velocity (the latter three 
depending on r and t). 

The formula can be simplified in practice. Emanuel [1998] took representative constant 
values for the density and the drag coefficient, p — 1 kg/m 3 , Cd = 0.002, and assumed a 
simple velocity profile, v(r, t) = v m (t) f (r / R(t)) , with v m (t) the maximum speed across the 
storm at time t, R(t) some characteristic radius of the storm, as the radius of maximum 
winds, and / a scaling function that is the same for all storms. Denoting the spatial 
integral as I = J 27rf 3 (u)udu, then 



Emanuel does not provide enough details about the integral J, but from the values of 
power, radius, and speed that he uses as an illustration, it has to be I ~ 2, and then 
k = 0.004 kg/m 3 . 

One can get an estimation of the dissipated kinetic energy E just integrating P(t) over 
time. Introducing an averaged radius over the storm, R m , and the PDI discretization 
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[Emanuel, 2005b], 

E = J P(t)dt ~ pC D IR 2 m J v 3 Jt)dt = kR 2 m PDI. 

The best-track records still do not provide the information required by this equation, as 
in general the storm radius is missing. However, for the Northwestern Pacific, since 2001, 
some radii are available [Joint Typhoon Warning Center, 2011]. Out of 11358 6-hour 
records for the period 2001-2010, 1671 provide non-zero values of the radius of maximum 
winds and the eye diameter. Averaging these values of the radii of maximum winds we 
get R m ~ 30 km. If we average the square of the radius of maximum winds and take the 
square root, then we get as a characteristic value R m ~ 35 km, which is the value we use 
for the North Atlantic, then, kR 2 m = 4.9 • 10 6 kg/m and 

E ~ 4.9 • 10 6 PDI, 

which yields dissipated energy in Joules if the PDI is in m 3 /s 2 . 

Under these approximations, as the North-Atlantic individual-storm PDI ranges be- 
tween 5 • 10 8 and 2 • 10 11 m 3 /s 2 , the estimated dissipated energy turns out to be roughly 
between 3 • 10 15 and 10 18 J, i.e., between 0.6 and 200 megatons (1 megaton = 4.18 ■ 10 15 
J). Nevertheless, the real range of variation will be larger, as the variability of the radius, 
which we have disregarded, increases the variability of the dissipated energy. 
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Figure 1. Annual number of tropical cyclones, hurricanes, and major hurricanes in 
the North Atlantic as a function of time. The data has been smoothed with a 1-2-1 
filter applied twice. Unsmoothed data are also shown for the number of tropical cyclones 
(tropical storms plus hurricanes). Corrections are not performed, neither in the counts 
nor in the intensities. For the original plot see Goldenberg et al. [2001] or Webster et al. 
[2005]. 
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Figure 2. Annual accumulated PDI (APDI) in the North Atlantic, both for the original 
series and for the smoothed one (by means of a 1-2-1 filter applied twice). No corrections 
have been applied. The mean of the original series is also shown. The year with the 
largest APDI is 1950, but followed very closely by 2005. The maximum of the smoothed 
series is for 2004. For the original plot see Emanuel [2005b]. 
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Figure 3. (a) Empirical probability density of 1966-2009 North- Atlantic individual 
tropical-cyclone PDI (494 events), together with a power-law fit between 2.2 • 10 9 and 
6.3- 10 10 m 3 /s 2 with an exponent r = 1.16. See Corral et al. [2010]. (b) Same data and fit 
as in the previous panel but for the (complement of the) cumulative distribution function. 
Notice that the power-law distribution is no longer a straight line in a log- log plot, due 
to the bending effect of the upper truncation, which adds a constant to the power law 
[Burroughs and Tebbens, 2001]. In this plot, the KS distance turns out to be proportional 
to the maximum difference between the empirical cumulative density and its fit (rather 
small). 
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0.43 



67±5 % 



e 1. 



Parameters 



23.3 

of the maximum likelihood estimation and the KS test for 



the PDI of the 494 North Atlantic tropical cyclones of the period 1966-2009, using the 
generalization of the method of Clauset et al. [2009] introduced by Corral et al. [2010]. 
The condition b/a > R m in is imposed. N ab refers to the number of tropical cyclones 
with PDI value between a and b; CN ab /N is the constant of the power law that fits the 
empirical distribution between a and b when the latter is normalized from to oo, its 
units are (m 3 /s 2 ) r_1 . The uncertainty of the exponent can be obtained from its standard 
deviation, a T [Aban et al, 2006], which, in the limits of large N ab and b ^> a, is given by 
o> ~ (r— 1) / y/N ab . d is the minimized KS distance, and d\/N ab a rescaling of that distance 
that should be independent on N ab , for comparison purposes, a and b are determined with 
a resolution of 30 points per order of magnitude. The p— values are calculated from 100 
Monte Carlo simulations. The last row corresponds to the period 1966-2007, analyzed by 
Corral et al. [2010], where an erratum was present in the value of CN ab /N. 
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Table 2. Same as the previous table but using the fitting procedure of Peters et al. 
[2010], which maximizes the log-range under the condition q > q m i n , where q denotes the 
p— value for fixed a and b. 
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Figure 4. Same North-Atlantic tropical-cyclone data than in the previous figure, but 
with a gamma fit for PDI > 2 • 10 10 m 3 /s 2 . The exponent is r = 1 + = 0.98 and 
the scale parameter c = 8.1 • 10 10 m 3 /s 2 . The fit seems reasonably good for 2 orders of 
magnitude. 
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Figure 5. Same North- Atlantic PDI data than for the two previous figures, but 
separated into two distributions, one for years with SST above the mean (large in the label) 
and another for below (low). The gamma fits, enforcing a constant exponent 1 + = 0.93 
(see text), show an increase in the scale parameter c around 50 % between cold and warm 
SST years. 



DRAFT 



October 5, 2012, 1:29am 



DRAFT 



X - 44 CORRAL AND TURIEL: VARIABILITY OF N. ATLANTIC HURRICANES 





SST 


MEI 


NAO 


AMO 


TlTC 


n H 


variable difference 
(PDI) high 
(PDI) low 
(PDI) ratio 
(PDI) difference 


0.49°C 

2.55±0.24 

1.79±0.17 

1.42±0.19 

0.76±0.30 


1.41 

2.09±0.23 
2.29±0.21 
0.91±0.13 
-0.20±0.31* 


0.79 

2.11±0.20 
2.26±0.22 
0.93±0.13 
-0.15±0.30* 
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1.45±0.45 


variable difference 

(DDI) very high 
very low 

(PDI) ratio 

(PDI) difference 
Table 3. Tor 


0.65°C 

2.70±0.29 

1.64±0.18 

f.65±0.26 

f.06±0.34 
) half: Mec 


1.92 

1.39±0.17 

2.10±0.24 

0.66±0.11 

-0.71±0.30 
m value of t 


1.15 

2.43±0.27 

2.22±0.30 

1.09±0.19 

0.21±0.40* 
re PDI for 1 


0.52 

2.66±0.31 

1.56±0.17 
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1.09±0.35 
iigh and loi 


9.47 

2.53±0.29 
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0.93±0.36 
n values of 


6.32 

4.26±0.41 

2.51±0.31 

1.70±0.26 

1.75±0.51 
seasonal 



averaged values of the SST, MEI, NAO, and AMO, together with the annual number utc 

of tropical cyclones, and the annual number of hurricanes, nn, for the North Atlantic in 

the period 1966-2009. For the latter case the averages are performed only for hurricanes 

(otherwise we could get artificial larger ratios due to the fact that years with high number 

of hurricanes and less tropical storms would yield higher values of the PDI). The difference 

of the variables between the mean for high and low years is shown, as well as the ratio and 

the difference between the mean values of the PDI. An asterisk denotes non-significant 

differences (at the 95 % level). Units of (PDI) and its differences are 10 10 m 3 /s 2 . The 

uncertainty of the mean PDI is the standard deviation of the PDI divided by the square 

root of the number of data. The relative uncertainty of the ratio is the square root of the 

sum of square of the relative errors of each value of the means. Bottom half: The same 

but for values of the SST and other variables above the mean plus one half of the standard 

deviation (very high) and below the mean minus one half of the standard deviation (very 
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Figure 6. (a) Same North-Atlantic PDI data separating this time by hurricane-season 
MEI above its mean plus 1/2 of the standard deviation (labeled very large) and below the 
mean minus 1/2 of the standard deviation (very low), (b) The same for AMO. (c) Same 
data separating by NAO > and NAO< 0. In this case the differences are clearly non 
significant. 
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Figure 7. (a) Same North-Atlantic data separating by annual number of tropical 
cyclones above and below its mean ± 1/2 of the standard deviation (labeled as very large 
and very low, respectively), (b) The same separating by number of hurricanes above or 
below its mean (large and low). This plot is not the one corresponding to the last column 
of table 3, as tropical storms are not removed here. Caution has to be present in the 
interpretation of this case, see text. 
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